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ABSTRACT 


An experimental and theoretical analysis was 
performed on low energy impact loading of thin (.940 
inches) graphite-epoxy plates. Six-inch sguar2 plates 
were subjected to dynanic impact. (below the ballistic 
range) and static loading. Th2 plate static strain 
energy and dynamic impact energy for failure *#vere 
equal and constant. Theoretical analysis was 
performed using both exact and finite element aathads. 
Small deflection th2ory was assumed and found to be 
inapplicable; the plates behaved in 2 Manner 
indicating that a lacjge deflection theoretical model 


would be more appropriate. 
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I. INTRODUCTION 


—- ee oe op ee ee oe 


Composite materials have a tramendous potential in the 
structural field. Although their future value in the 
aircraft industry is exceptional, there still cemain many 
unknowns. The response of graphite-2poxy composite plates 


to low energy impact loading 1s one of these unknowns. 


Low energy impacts are a common occurcrence on 
Operational aircraft. A mechanic dropping a wrench, or a 
stone thrown by a tir2 are typical examples. Ee sis 


important to know how composite panels react +9 these 
impacts in order to be able to confidently predict a damage 
threshold. An accurat2 ndieal for damage would be helpful in 
both the operational and dasign fields. Operatioaally, the 
ability to ‘assess damag2 to a panel which is known to have 
been hit, but which has ad visual manifestations 9f damage, 
will produce better maintained and safer aircraft. ne 
designer will also be better able to produce ai aircraft 


capable of operation throughout its design environment. 


A series of tests, both static and dynizic, were 
performed to analyze th2 response of a Simply supported 
plate to impact loadjng. The plate was a synmetrically 
laid-up graphite-2poxy plate six inches sqguars. fhe dynanic 
tests used a variable-+~mass steel p2netrator of spherical 
Shape, at impact speeds of ten, fifteen, and twenty feet per 
second. Each plate was inpacted by incrementally increasing 


masses until failure occurred. 


There were three basis types of failure expected. These 


included: static shear-sut, where tha penetrator sheara2d 
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through the panel causing only local damage; flexural, wheres 
the failure was causei primarily by bending 311141 membrane 
stresses; and baliistic, where complete penetration 
occurred. It was assunedid that wave stresses would have 
little influence, and that flexural Failire would 
predominate. Small deflactions ware initiaily assumed and 


membrane stresses were neglected. 


The KP EChort Hypothesis for plates was adopted: 
straight lines normal to an und2formed middle sucface were 
assumed to remain straight and normal to the nidile surface 
in the deformed state. Further, it was assumed that rormals 


to the middle surface do adt change length. 


The investigation proceeded in the following manner: 
first, the constitutive properties of the plate material 
were determined both analytically and experimentally, then 
the experimental apparatus was set up and calibrated to give 
the desired range of impact velocities. A number of 
non-destructive technigu2s for inspection or danaged flates 
were considered and tested. Pinally, a series 9f static ant 
lmpact testS was run, ani the rasulting damage was evaluated 
uSing both nondestructive and destructive test methods, 


including scanning-electron-microscopy. 


Concurrently a theoretical investigation was undertaker 
to see whether simple linear theory could be of use in the 
analysis of the plate impact problen. Static and dynamic 
solutions were sought using SAP IV, a E£inite-element 
structural analysis progran. Initially, known closed forn 
solutions were compared with SAP [VY solutions for static and 
dynamic loadings. fhen specimen constitutive prooerties and 
experimental loadirg paraaeters were entered in ti2 progran, 


and solutions were compar2i with experimental data. 


The procedur2s used and the results gained ace detailed 


is 





Himere  ZOllowing Sections. 
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ieee cINSOTITUULEVE PROPER Mises 


fav ArtRtAL Se Roe ERTIES 


The plates used were constructed From 12 inch wide rolls 
of UCC Thornel 300 (uatwisted) Sraphite Fib2rs, in a 
Rigidite 5208 matrix, sapplied as a B-stage preodreg by the 


NARMCO Materials Division of the Celanese Corporation. 


Twelve-inch Sjuare oan2els ware constructed following the 
Manufacturers recommended procedures, and thos2 of Linnander 
(Rer. 1]. These proceilures were designed to jive a 554 
fiber content by volume. The lay-up was eight lanina in 43 
Symmetrical configuration as indicated by Fig. 1. fhe 
Standard laminate orientation code for this lay-up [Ref. 2], 


was (0,+45,0) . After the proper cure, the oanels were 
Ss 


jJuartered to give Six-inch sguare plates. 


Some of the lamina properti2s predictei by NARACI 


Materials includeji, in pounds per square inch: 


Flexural modulus (longituiinal)=22.0 E6 

Flexural modulus(transvecse) = 1.6 ES 

Tensile modulus of coupon sample (Logitudinal) =20.5 E6 
Tensile modulus for coupoa sample (transverse) =1.47 E6 
Plexural ultimate stress (longitudinal) =295 E3 
Flexural ultimate stress (transverse) = 11.5 E3 


Tensile ultimate stress for coupon sample 





LAMINA FIBER 
ORIENTATION 





Xx 

LAYER 9 
. 

L040 in “45 
' 4 0 
CLOG abel 5 0 
ee el: 

Figure 1 —- OR TENTATIIJN OF LAMINA ON A LAMINATED PLATEe 





(Longiscwa@2na lL) —2 15 53 

Tensile ultimate stress for coupon sample 

(transverss) =6./5 E3 

Compression ultinate stress (longitudinal)= 198 8&3 
Interlaminar shear ultimate= 17.8 E3 

The published properties aid not include either of the 
Poisson ratios 2r the sheac modulus. Th2 nominal thickness 


of a laid-up plat2 was .J49 inches. 
See lemon oeen ome ows FErULEVE PROPERIIES 


Jones [Ref. 3] developed equations for the constitutive 
properties of a lLlaminat2 from known Lamina properties and 


laminate lay-up. 


Beginning with the known properties the constitutive 
properties of eaca lamina in relation to the laninate w#as 
determined. Lamina with Fiber direction coincidjdeat with the 


axils of the plate; and plane stress: 


és = Q E 


KIN, AYER K K 


For laminas that were not on the principle plats? axis the 9 
matrix was Todt ed With I Suitable coordinate 


transformation. [his applied to the +t 45 degree Laminas: 


Q, = QuCOS'O + 2(Qi2 + 2 Q,d\SINOCOSO + QnSIN'6 
Qn (Op + Qi2- 1Qee) SiN*6C05B8+ Q2(SIN'O +0819) 

Qn = QuSin"O + 2(Qi+ 2Qcc) SIN'ECOS*O +Qx2C0S 6 

Big = (Qu-Qi- 2Qce) SINGCOSS + (Qy2-Qz2+ 2046) SINOCOSE 
Dae = (Qu -Qir -2Qe0) SING COSO+ (Qiz-Qir+2.Qee) SIND COS'S 


ey 





Joe = (Qi +Q2 -2Qi2-2Qce)SIN‘BCOS'O + Qe (SINS +COS1O) 


Foc small strains and linear elasticity, utilizing 
classical laminate theory, the stresses in any layer were 


expressed in terms Of midile surface strains ani survatures. 


fol Jetted +206 5, 


By integrating the stress in each lamina over th2 thickness 
of the laminate: 
Ye 
N O, dz 
X “Ib * 


and 


Y 
M, )_ eX dz 


the following relations w2re optained: 


in} = [a] te} + Le] 
wg = [oe] fef +o | 


where for extensional stiffness: 
n 

=> (Q. Z ~-Z_) 

K=| (Q;, ‘“ K Kl 


BOG CompimmiG Stafinass: 


2 2 
B =, 22 (ZZ) 


for bending stiffness: 


= yt 
H 


{ 


The plates as constructed ware synmetrical. fhis forced 


ho ’By  CbeCoOmD lida t © ins tO 22050. 


To facilitate calculation it was desirable +2 express 
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EhemmcovsertUtlVe properties in “the @orm of snginesring 
constants of an octhotropic material. The finite element 


program used woild accept data in this form only. 


The Design Manual (Ref. 2] supplies a means to 
approximate the OUtQoeL ODI Cc sonstitutive properties 
utilizing the A, or extensional stiffness, matrix; 

2 : 
E xx = (Ay; Az2-Aj2)/ (THICKNESS)Aoo 
E yy =Exx Aza Ay 
Gyy = Aa3/ THICKNESS 
V yx *Al2 ZAn 
Vxy 7 A227Aj2 

A Simple computer program was d2signed to take both the 
lamina constitutive propacties and Laninate constcuction as 
Poput.. The program, in BaSic Language, outputs the 
following: the indiviiual lamina stress-strain properties 
in celation tos tha lLlaninate, the 2xtensional natrix, the 
bending matrix, the orthotropic constitutive properties, and 
the orthotropic stress-strain relationship. This program is 


Appendix A. 
Soe EAP ERP RENTAL CONSTITUTIVE PROPERITES 


To verify the calculated constitutive properties, a 


series of tensil2 tests was performad. 


Three? specimens wees preparei in accordance with the 
d2sign manual [Ref. 2]. [aey were nine inches long and one 
inch wide. The first Specimen was cut from a twelve-inch 
panel with the zero degjyr22 fibers ranniny parall2l to the 
long dimension. The s2cond was cut so the z23ro degree 
fibers ran at a right angl2 to the long dimension. The last 

S 


specimen was cut so one of the forty-five degr22 lanina' 
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fibers were parallel to tne long dimension, and the zero 
degree fibers were forty-five degr2es off the Long specimen 
axis. One-and-a-half-incn long and one-inch wiie shoulders 
were glued to both sides of the outer osne-and-a-1alf inches 
of the specimen in order to distribute the local stresses 
caused by mounting in th2 test machine. The eiges, where 
the shoulders met the spesimen, wer2 beveled to l2ss than 30 


degrees. 


BBH Electronics SR-4 Strain gage cosettes were installed 
9n each specimen near th2 center. Thes2 rosett2s were type 
WeekR=zon-lesGo Witn 120 San resistance and 2.04 jgaje factor. 
They were wired intd a standard Wheatstone bridge strain 


measuring circuit. 


The three specimens were azdunted in a RIEALE testing- 
machine. As they were slowly loaisd to failuce, strain 
meaSsurementS wer2 inccenentally takan. Figures 2, 3, and 4 


are plots of the test data. 


The orthotropic SXRSt CUCL Y= peodperties can be 
determined from the graphs. From this data, the Design 
Manual {Ref. 2) gives a simple calcsalation pcodcedure to 
datermine the shear modalis: 

U=(E, + Eqxt2Vi7Ey)/ 8Ul-VipV rn 
5-2 ¥ Ex 
LT SU-Ey, 





From the graphs: 


Fis 11.65 E6 
Er= 3.1 £6 ‘Nir 52 
Ey= 61.6 E6 Vru*.-l4 
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PO vVeLlif£y Poissons ratid: 


E 
a (a ) = 138 % .14 


The shear modulus calculation gave: 


S122 2OG jee 

Summary of constitutive properties from tests: 
Ei ll.6 E6 
E-= 3.1 E6 Vi 77-82 


Gis= 2.1 E6 Vrit4 


Using this data the str2ss-strain relationships will be: 


12.5 1.75 O 

O? = .75 334 0 E6 E 
O O 2.1 
i= 


PeeeGOnePAnRTSON OF CONSTITUTIVE PROPERIIES 


Since Poisson's ratios and the shear nojulus for an 
individual lamina were 19t given, these properties were 
estimated and entered into the constitutive properties 
SonpUber program. Using the orthotropic properti2s obtains 
Prom the tests, the program was used to iterate for tne 
unknowm lamina properties. The results ware disappointing. 
The calculated aodulus i2 the ¥ direction was eight percent 
Pugteonesnedar Modulus #is forty percent high, and the 
Poisson's ratios were thirty five percent high when the X 
direction calculated modulus was forsed td equal the test 
data. 


AVSUrve eOtmememoesign Manual (Ref.2] indicat2i that tue 
constitutive properties obtained experimentally were 
comparable to conmon, medium strength graphite epoxy places 
Of SSinilarapcomceric ti Dh. The orthotropic values obtained 


from the test data were ised unless otherwise note2d. 
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E. FIBER-EPOXY VILUME FRAZIION 


Two independent tests were conducted to verify the fiber 
and epoxy volume fractions. Ine method was a chemical 


technique, the other was a quantitative nicroscosic method. 


Handley and Cross ;R3af.4j sutliaed the chenical method 


used. Two fritted crucibles and two samples were dried in 
an oven at 180 dejrees c2ntigraie for thirty ninutes. The 
crucible and coupon w#eare then weighed. [Th2 coupon was 


placed in 80-90 degree centigrade, concentrated nitric acid 
EOL fifteen to twenty minutes. The hot ficriecmpacid 
dissolved or "digested" tae epoxy matrix. Aftar the epoxy 
was dissolved, the cenaiaing fibec and acid solution was 
rinsed with a vacaum assist in water and acetone. The 
crucible with the fibers was dried in an oven at 120 degre2es 
centigrade for on2 to two hours, ani then weighed. The 


results were: 


Sample One Sanple Pwo 
eon 2 specimen weight ma gcams 3525 
14.1733 specimen and crucible (before) 14.2385 
13s 9225 specimen and cruciple (after) Aseovo7 
~2508 epoxy w2ight in graas eo was 
23.6% epoxy by weight 24.8% 
e255 specific Jansity ait 2poxy ey Bers. 
Ug Se Speciiic Jensity of graphite es: 
94-68% fiber by volume 63-67% 


This compares favorably with the desired 65%. 


A quantitative microscopic analysis was done t> further 


verify the fiber/fepoxy mixture. 


A section was taken from a plat2 using a diamond cut-oref 


ZS 





Saw, aha Cold mounted in plastic. [t was than polisned 
isngegernaing belts of 5) and 320’ grit, smory piper of 9 
and 3/0 grit, and polishing cloths with aluninumn-oxide 
Slurries of 15 and .05 miscometer grit. The misrostructure 
showed sufficient resolution with just polishing, so etching 
Was not used. A nitric atid etch wouli have selectively 
attacked the apodxy, if iesiced. The sold mounting process 
caused slight sample degradation. Th2 plastic attacked the 
epoxy in the outermost layer. The attack was very slight 
and confined to dissolviajy the extrame surfase layer of 
epoxy. If there was a well davelopei crack, the plastic 
would also f111 these with a corresponding light attack on 
the surrounding matrix naterial. The plastics mountioag 
material was also Significtantly softer than th2 graphite 
epoxy plate making it iifficult to obtain square edges on 


the sample. 


Figure 5 details the ceasults of the gqiantitative 
dJialysiSwenem tla StEratiog 15 EhevsSIX nicrophotograph used 
for tire AMAL ES WALeb 2S a4 picture of the middle two 
layers of a plate. The i32ta indicates a fractioaal length 
Of linewWethineestme  f1DSt aS .597. This equates tos the 
volume fraction giving a 55.7% fibec volume, whith concurs 


with previous analysis. 


The average dimension of the fibers was also determined 
froma the mMicrophotogcaph USi0deeeseemndard quantitative 
microscopy tecaniques. The average dimension was 5.4 
microns. Since the specinen was cut at a right aagle to the 
principle fiber orientatioa, this dimension gives a close 


approximation of the average fiber diameter. 
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Tit. EXPERIMENTAL APPRATUS AND P2SLIMINARY TESTS 


A. EXPERIMENTAL SET-UP 


1. Penetrator 


Two variaole mass p2netrators were constcucted and 
ace illustrated in Fig. 5. The prinary penetrator was the 
lighter, slimmer one. [t was sconstructedi with a plastic 
pipe body sized to fit in the launcher tube. The mass could 
be adjustei from approxin2aatly .014 to .993 slugs. This was 
accomplished by varyiny the weight from seven ounces tod 
three pounds with lead shot. The Leading podinee-of s ehe 
penetrator is half-spherical shaped and made of steel. The 
larger all steel penetrator had a mass which coulji be varied 
SEaeting at ».0oo SLUGS." Tals penetcatoc had the capacity to 
change both the leading point's shape ani its construction 
Material. The second o92n0etrator was not naeeded in the 


dynamic tests. 


A MasSive steel fixture was built for both the 
Static and dynamic tests. This fixture is illustrated in 
Fig. 7. It was carefully machined to have a flat test 
surface and to have th2 top and botton surfaces parallel. 
The fixture was made fron heavy steel channels with a 


One—-HhalE=1nch toiecwmestee!, base, 
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Figuce 7 - PLATE TESTING PIXTURE 
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The square test sudports measurei four-and-one-half 
inches on the insid2. Using a six-inch plats, 2ach plate 
edge over-lapped the supd drt by thr2e quarters of an inch. 
The fixture also hai an additional frame and a series of 
clamps to secure all edges of the plates to the fixtures; 
this frame would prevent all edge rotation. Plat2s could be 


mounted and tested in either th2 simply supported or ciamped 


edge modes. Ihneve@e | Veber Tiere fale sls cal jie) eee: epbGlists) epg ie lets 
fixture to prevent any pressure build-ups duriagy dynamic 
displacements. 


Ss Launcher 


The desirei impact velocity was achievai by simply 
dropping the penetrator from an appropriate h2ight. size) 
insure that the penetrator would hit the desired spot on the 
plate, and to provide for cepeatable, easily aijustable drop 
heights, the launcher sadwn in Fig. 3 was constricted. The 
launcher could be varied in height above the plate via 43 
screw crank for fine velocity adjustments. Both d2netrators 
were just slightly small2c in diameter than the inside 
dimension of the launshec tube, which was usel td) stabilize 
the penetrator's flight. Launch was initiated by a switch 
Bontroalled solenoid @aich cenovaei a bar holding ths 


penetrator in place. 


A frame was constructed out If steel scarfolding, 

I1C ei Ustrateds in Fig. 9. The frame was plumbed 

VCCCICa wancmconcermimcted 22 a heavy-load cated Concrete 

rel oo an The launcher assemply was moved up ani down the 
3) 


frame by a pulley and coo2 arcangenent to Obtain course 
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velocity adjustme2nts. 


32 








Fanci: e. 3 


PENGERATOR LAUNCHER 


33 


hy 
Ez 
af 
ho 











Figure J > -LAUNSHER FRAnS 


34 





5. Static Tests 


Static tests J were conducts] on a RIEALE Cectung 
machine uSing the base fixtures, plate and appropriate 
penetrator. Deflections were measured Wil tUh@eaes dial 
displacement iniicator fron the cross acn of the testing 


fiaenine to the working surface. 


Be. CALIBRATION 


The launcher-£rame coabination was calibratei for both 
impactors at three different speeds. Ten feet per second 
Was the lowest practical speed obtainable due to Launch tube 
interference. Twenty f22t per second was tne highest speed 
desired in order to keeos in the flexural failure nmoie. 


Fifteen feet per second was also used. 


VremcaelbratiLem Was’ done electronically. We isles 
Sensitive photo-transistors and their power supplies were 
placed two inches apart, bracketing tare plate fixture's 
position relative to the launcher. Iwo ligat beams were 
focused across the penatrator flight path ones the 
transistors, and an elestronic timer was placei setween the 
transistors. The wooden fram2 supporting this aposaratus was 
Painted flat black to prevent spurious imputs. Speeds 3£ 
ten, fifteen, and twenty Feet per second were salibrated 


Ret hin Gewe percent. 


oye, 





eee ACT DURATION 


pieowrime, ands Strain-adpilituds aistory of a typical 
impact was obtained t) assist in the theoretical 
calculations. A graphite epoxy plate similar to the teste] 
plates was fitted with a strain gage near the center of a 
face. The strain gage was connectei to a Wheatstone bridge 
and power Supply. An ossilloscope was connected across the 
bridge and set to triggec on impact. Ine scope was set for 
Single Sweep moje, ani had a camera mounted tos record the 


traces. 


Fi gue NC eset) SCOpe picture taken using thes larger 
impactor at its lightest weight. The speed of impact was 
approximately ten feet per second. fhe oscilloscope time 
scale for this photograph was .5 miliseconds per tentimeter. 
The amplitude waS not calibrated and gives relative 
information only. The stcain in th2 plate reached its peak 
ValucweeOneuwmibliS2cOldueaeccOrding t> this photograoh. 
After the impact, examination of the plate deternined that 
the strain gajge leads had separated fron the gage leaving 
the validity of the results in doubt. The test did give 
several useful pi2ces oF information. [It gave an indication 
of the lmpact's initial saape and ris2 time,and, since wave 
Stresses would have reataed maximum values in nich quicker 
time ranges (microseconis), the test confirmed the 
assumption that a flexural or menbrane mechanism would be 


the primary loading mode. 
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Pree vou-vBhs TRUCTIVE TESPENs 


It waS necessary tS exaMine each plate after every 
PpterelOading to deteravlis Mf any dazage had occirred. ihe 
tests had to be guick, in2xpensive and non-destrustive since 
many tests would be ceqyuired during the course of the 
experiments. Several nethods of nadn-destructive testing 
were considered, including: ultrasonic sound, X-cays, and 3 


new liquid crystal process. 


APA womselone Lidustrisas “[as., Sperry Division, 
Pic ue iMenerlectoscope Ultrasonic Pesting %1cChminep vas 
examined for its ability t>5 meet the requirenents. The 
machine was eguipped with 2.25 MHz ani 5.0 Mdz zero degrees 
maanrsducers., 8 ft also had’ 2.25 “Hz transducers with sound 
Ai Ut angles of 45, 50, and 90 degrees. All the 
transducers, including vacious combinations, wer used 0a 
two different graphit2 2poxy panels with Known EClaws. Ona 
panel had the standard sucface finish, and tn2 other was 
Specifically manufactur2i tS give a very snooth surfaces 
finish. Both the shear ani through anodes wera attenpted, 


along with the different angle input modes. 


The panels with nocaal surfac2 roughness, fron th2 
bleeder plies used in nanufacture, ware difficult to test. 
[he S@@memereturn £Erom both the front (input), and b32ck 


SULEACeomMtomPost 1h tas sluttec caadsed by the coighness. 


The plate with a snooth surface allowed a reasonable 


input and back surface echo, but the plate was s> thin that 
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these rceaeflections occupied the entire useful canjye of the 
echo returns. Knoown Flaws could be Jetasted in the clutter 
uSing the highest frequency, but it was J2ifficult to get any 
continuous, consistent ceturns. This was tru2 even with 


known, large flaws. 


Other modes of th2 machine were tested with negative 


results. 


Fonda Dlhateweas “thin weas F11s;asnighss re oquene, 


transducer might have given useful results, but for the 


purposes of this inv2astigation ultrasonic t2sting was 
discarded. 
2. X-Ray Testing 


A portable X-Ray Wachine was triei as a means of 
non-destructive testing. A General Electric LX¥-40 Portable 


Industrial X-Ray Jnit was set-up to examine a gcaphite 2apoxy 


panel with known flaws. Polaroid Eilm was exposed at 
Various kilovolt and tine settings to daternine the 
procedure's usefulness. Even at the Lowest setting of 70 


Kilovolts and 10 seconds, tae plate did not absorb enough 
energy to produce an inage. The pictures obtained had no 
contrast and no detail. A nachine capable of evan smaller 
ejergy output would probaosly still have difficulty obtaining 
emyecontrast due to the light atomic weight of the carbon. 


This means of non-destructive test was abandoned. 


Seystal fsstang 


Shaun { Ref. 5} develop2i a methoi for using 
temperature sensitive ligqiid crystais for damage evaluation 


in gctaphite epoxy laminates. This mathod was used with 
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excellent results. 


After Oreparing plates for testing, they were 
painted flat black on 922 side. baud crystals, wiaeeh 
temperature sensitivity of 30+1 degrees centijyrade, were 
mixed with water and a snall guantity of wettinjy compouni 
(liguid detergent). This mixture was then Spray painted on 
the plate to be tested. 


After each impact during the dynamic t2sting, the 
plate was examined for danage. This was accomolished by 
selectively heating all, or just a portion, of the panel, 
and watching the heat traaisfer cause the Liquid crystals to 
change color. Any damaged area would experience a different 
neat transfer rate ani thus would change color at 4a 
different time from an andanagei part. This method was 
sensitive enough ED show 2ven Slight sompressive 
indentations in the sucface prior td. their becoming visually 


or tactually evident. 
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iV. Sc iNi ie SLCSUENT SANALYSES 


A. FINITG ELEMENT MODEL 


Sane EVye a —“Strucrucal™ analysis program for Stateerand 
dynamic response of lineac systems oy Bathe, Wilson and 


Peterson [Ref. 5] was used as an analytical tool. 


Two finite elament adsials of a plate were made for 
analysis. The first model was a coarse nodel coasisting of 
only 16 square plate elenants, 1.125 Lathes on a side, for a 
meme 2och saiiare plate. This model was used for initial 
Mcoarse™ investigation at a small cost ih computer tim:. 
The second model was a fine nodel with 324 piat2 elements, 


-25 inches square, and us2i considerable computer time. 


The results produced by SAP IV using these nodels were 
compared with known, exast soiutions foc aluminiua plates. 
Timoshenko and Woinowsky-Krieger [kef. 7} solved the simply 
supported, concentrated load problen, which was used as the 
rererence SOLUEL On. Piguce 11 shows a plat “etter 
d2flection versus position 9n a plate (from an 21328 to the 
center keeping Y sonstanot) under a 109 pound 15ad. The fine 
finite element model is .4% stiffer than the exact solution, 
Where the coars® moiel is 95.5% stiffer thaa the exact 


SO2utlLon. 
The dynamic analysis section of SAP was addressed tod 


compare its modal freygiency with an Xia ole SO Wate Oke 


Meirovitch fRef. 3] gives an exact (iso luc o ieee o eco 


4] 





vibration frequencies of a rectangulac plate: 


2 2 29 
F REQUENCY = Vo (2Y¢(PY myn = 1,2... 


iitemecoatse E£ilnite elaneat model results were within .-5% of 
the exact solution. Sins2 this agreement was considered 


Wome me ial satisfactory, tas fine molei was mot cua. 
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Bemeot nk kooe OD LSITR LEULLON 


Dime onder tO Ttmit the required anoint of conputer d2t> 
reduction by determining where initial failure could be 
expected, an analysis was perforned using 3AP IV to 
determine the stress jistribution on a typical plate. 
Figure 12 illustrates th2 results of this study on a sinply 
Supported plate with a concentrated Load. The stress 
concentration w#as the highest at the center of the plate, 
and this indicatei that initial failure analysis couid be 
limited to the immediate vicinity of the caiater of the 


plate. 
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ce FAILURE CRITERION 


Tests made by McQuillen, Llorens, ani Gause “2af. 9] on 
composite beams demonstrated a flexural failure aodde. Thus, 
For this investigation th2 assumption was maidie that the 
flexural failuce mode woild predominat2. Since 2ach lamina 
waS Significantly weaker in the transverse direction than in 
the direction parallel to the fibers, and since the 
Sutermost layer would hav2 the highest stress l2v2l, it was 
expected that the prelininary failur2 would be in the 
Qutermost layer i1 the transverse direction. After this 
failure, the load normally carri2i woald be traasferred to 
Mmeemgext Layer. Lf this Lanina could adt carry tne Load ain 
Pe met hancoveobscmaibccti1on, Eaillura 3f£ ths plat? wis accancd. 
Bach lamina's longitudinal strength was so larg2 that a9 


fagture Was allowsd in that direction. 
Retina tOenones I22f. 3], the Strains are: 


4 7 Ser aaa! Q Uo 
Ey: $e z os VW SER 


of 








: Ji, 
¥-EL AEN 


Since it 1s assum2d that ta2 length sf th2 midile surface of 


the plate is constant, tha first partial terns can o0De 


oy guy = fet 


the middle surface curvatures are found. Having these, and 


dropped. 


Using, 


Sueommcrariness matrix fot each lamina, the lania2a stresses 


are found fron: 
o| -| @ | € E “24 
K K 
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Knowing the stresses ani the ultimate failure stress, 


failure can be predicted. 


A computer progran waS written in Basis Language to 
evaluate failur2. The insut consistei of the moments per 
Guee Length ina the X, ¥, and XY directions. [h2S prograsn 
Started by computing the stresses in the outer two Layers, 
as outlined above, using appropriate coordinate 
transformations. it then set the transvers2 and sha2ar 
modulus for th2 outer Lamina t>2 Zzerd and computed the 
stresses on the n2axt Llay2r. Failur® was assum3i tO occur 
when the transverse ultinate failure stress was axceeded in 
both the outer and next layer. The program is iasluded as 


Appendix 3B. 
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VeneoOLALEC TILES! RESULTS 


Seven statics tests wece performed, six involving loading 
wom carlure. Three tests used the Lacger penatrator and 
Simply supported boundacies. Two used eves snaller 
penetrator and were also simply sipported. J1e test was 
conducted using fixed boundary conditions ani the small 
penetrator. The laste plate was not Loaded until Faun. 
instead, it was partially loadei and then used LOG 


microscopic exanination. 


Figures 13 and 14 show the cesults of all static tests. 
The static failure energy, as outlined by AYAcQuilien, 
Llorens, and Gause “Ref. 39] for beams, was one half the load 


times the displacement. 


All failures, with on2 exseption, involved the same 
failure mode, and Figure 15 and 15 are typical 2xamples of 
this mode. The only eaxtaption was the test with fixed 


sdges, and Fig. 17 illustcates its aoile. 


Mimengeche Course of the Static t2stS it Was “Geman 
that a "yield" point was expecient2d. fhis p9iit occurred 
at approximately half the ultimate load when the panel wouli 
experience a decreas in” load with no 2a¢rea ss era 
displacement at the center of the plate. Pigure 14 
mraecates that this was comnon td both peaekrators and 
Dolg@dany conditions. <Ionsequently, a plate was Loaded to 4a 
Popememmst beyond this "yield" point, and then ssstrened and 
prepared for microscopic 2xamination. AS th2 plate was 
loaded a large deflection occurrei in the plat2 which was 


Very Similar in size and Shape td the stardari failures. 
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Figure 18 indicates this ieflection area, and shows were the 


sections (marked 5) were taken. 
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Piguce ) 18 Pear oNs 3 SRCTIONS REMOVED GIR REECROSCIPTe 
EXAMINATION 
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Vi. JINAMESTIEST ntovukes 


Figure 19 simmarizes the results of the dynamic tests. 
Each plate was impacted c2apeatedly by the penetrator, at 
constant velocity but with progressively insr2asing mass 
until failure occurred. Failuc2 was evidenced ovpy a loud 
cracking noise, victually no penetratoc rebound, and gross 
delamination. The energy c2qguiced for failure is indicated, 
and 1s one-half the mass times the velocity sguaced. After 
each impact the plates wer2 exanined for damage using the 
liquid crystal iaspection technigue. In all of these tests 


the smaller siz2 impactor was used. 


Each failure, with one 2xception, showed the same 
Failure mode as the static case shown in Fig. 15 and 16. 
The one exception occtucred at 29 feat per second, and 
approached ballistic in nature. Figuce 29 shows this 
failure with its almosc complete penetration, aad splaying 


of the delaminated layers on th2 back (non-impact) side. 


Early in each test, damag2 was observed in the plates. 
This damage was lLocaliz2i in the immediate area of Enea st. 
When it started it was oarely discecnible using the ligaid 
ceoystal technique, but as the tests projressei it bdecane 
easier to see. Aft2r nanny impacts, the danagead area was 
dJ2laminated, and this coild be determined by careful tactile 
examinations. It did not appeac to have any effect on the 
the ultimate failure e2nergy. This "bruisei" area was 
subjected to a microscooics exanination. It was section 1aAe 


in Fig. 18 and is shown enlarged in subsequent figures. 
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Vil. DESC USS LINE ea eo UL ho 


There were three striking phenonena encountered during 
the course of the tests. [he first was that the str2in 
energy for railure in tha statics tests and the imoact energy 
for failure in the dynamic tests were egual andi constant. 
The second was the initial dynanic ianage that occurred well 
before the ultimate failure. The third was the large 


deflections encountered. 


As the tests and analys2?s progressed, the large 
deflections had an important influence on the sorcrelation of 
experimental and theoratisal results. The plates were .040 
Mmenes thick ani statis failure oscurred at plate centsr 
teplections Of  approxinately .35 inshes,. cis means 
deflections were eight to nine times larger than the 


thickness of the plates. 


An analysis was attemptei to correlate the actual 
results with theoretical solutions. five exact and finite 
element solutions were dJotained for th2 statics case and an 
intensive structural dynanic response analysis was attemptei 
for the dynamic cases. Thes? investigations wece thwarted 
by inadequacies in the basic assumption oye small 


deflections. 


The plates all 2xperi2nced a comnon failure node for 
both the static and dynanic tests. YWicroscopy was utilized 
to examine this ultimate failure sairface and aechanisn. 
Furtherc areas addressed with microscopy included: a 
pre-failure dynamic danag2 zone, or "bruise," that developed 


Gnmeecach “panel: 4 "yisldi® polnt thal occurred aueing staers 





ls5ading at approximately 222 half the ultimate load; and the 
fiber failure mechanism. 


A. STRAIN RATE DEPENDENCE 


AS indicated in Fig. 14, in ona lLarge-penetcator test 
the specimen se2med t)3 od2 significantly stiffec than the 
other plates, all Ot which Showed quit2 Similar 
load-deflection curves, cegardless of panatrator size. The 
only parameter which diffec2d in this "stiffer" t2st was the 
crate of loading, which was highec than that 1sed in thea 
other tests. No furth2r gualitativ2 analysis was performed, 
but the rate of loadiay was kept as low as possibl:. 


Further tests thea yieli2i consistent, repeatabl2 results. 
De FAILURE MODzS 


All failures, both static and dynamic, showed the same 
fFailurc2 modes (Fig. 15, and 16]. There were only two 
exceptions. The first was during testing of tha fixed plate 
with the smaller senetratoc. This failur2? was illustrated 
Pie Fig. 17, and was 3 forn Of Snhaar. AlthSugh ths plawewyas 
stiffer than the simply supported plates, it responded in a 
Similac manner until th2 "yield" point was reacaed. After 
that point, the panetratoc virtually sheared through the 
plate. The second 2xr2ption was on2 test at 20 feet der 
second. The failure mod2 approached ballistic in nature. 
Figure 20 illustrates the almost couplete panatcation, and 


the splayei lamina arouni the exit 4aole. 


The -choice of 20 f2Sst per second as ayia. fone 
t 


Gymanic tests was JusStiftisid. Since ofe of ths 2 ue cs 
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conducted experienced this failure moje, this speed range 
Hesein the transition £con flaxucal failures wwede =o 4 


ballistic failure mode. 


The size of the plat2 was sufficient to keap the plate 
in the flexural failurs, insteai of shear failure, mode 
during static tests. All tests in the simply supported noise 


were flexural in nature2. 
1. Uhleroscopre = uxalitatuonuss 


Section” 1A2 [Pig. 13] was cal andepounteimser Light 
microscope examination. figure 21 shows the ragion outsiie 
of the gross delamination area typical of failure. This 
section appears relatively 1 Pest with vary Minor 
delamination. [There ace2 several small lamina fractures in 
the outer two lay2rs cunniog with the fiber dirastions in 
the matrix, and thes2 ninor cracks ar2 associated with 2 


slight delamination. 


Figures 22, 25, ini 24 show the Lamina i1 the gross 
delamination regidn. Th2s2 show extensive delanination, and 
numerous large cracks if the Watrix, between fibers, 
penetratinjy entire Jlamiaas. Bach dealaminatioa betw2e2n 
laminas appears to have at Least one crack penetrating the 
lamina directly above or o2alow. Figure 24 illustrates this, 
With delaminations ani issociated fractures of tae adjacent 


laminas. 


Totally flexural failure damage shoulji have been 
concentrated in the out2acmost layers where tae Largest 
Stresses occurred. Figjgucas 22, 23, and 24 ali denonstrate 2 


Srgniticant failure of the middle o tye layers Ripe oe 


61 








Saw 


aN ae 


< 


~ * 
AS 
ae 





Figure 21 - LISHT MBERISIIPE, FALTLURE SURBASS eee ole 
SGT BON J Aas peel Joe 
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Figuee 22 - LISHE MICROSCOPE, FAILJRE SURPACE EX SORE, 
SECT RON Ase Ux 
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Pigure 23 - LIisHr MICORISCOPE, FALLURE SURES s SXKAMINATION, 
Se SEEON 1 aa ox 
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Figure 24 —~ LISHT MITRISCOPE, FALLURE SURFACE. s Giessen 
SECTION@AAS 2.2206 








transverse di Gee einen The assunption that the middle 
Surface of the plate iid not 2xperience a2xtension, 
therefore, is invalid. TInis would iniicate that theoretical 


analysis based on small daflection theory is inappropriate. 


In Figuces 21, 22, and 24 the danage is nost intense 
on only one side of the plate. The outer two lay2rs on the 
Side opposite the applied load appear almost intact. A 
quantitative examination 2f all failed plates showed an even 
distribution of damagei outer layers. Some plates showed 
gross lamina fracture in 39th outer Layers. Jtaer panels 
had damage on only the Loaded sii and vice-versa. Thus, 


tiewtwo Intact layers acr2= not considered significant. 


A section was cut for th2 scanninzy electron 
microscope (SEM) evaluatioa as indicated in Fig. 13. SLoice 
this waS an area of gross delamination, the top (or loaded 
side) two layers separatei during sectioning. fhese Layers 


were prepared for evaluation. 


Heowie 2o shOws A typical, fess tus, won wel ac wee 
Spmexy) matrix cunhing tacodugh ths Sutee Layer: Figure 25 
Jeamonstrates a sinilar fracture, further illustrating the 
complete penetration of tae lamina. fh2 antire suter layer, 
a zero degree layer, was extensively cracked. FIZUreS Weer, 
28, 29, and 30 show typisal fractuc2s in the first 45 degrees 


layer. 


Figures 27, 29, and?30Phave cGemnants Of Stroman 
inner layer (45 jiegrees) which are evident as spoxy and some 
fiber fragments running at fight angles t> the original 
layer. [his was typical of ally the ) gross (istaniqieren 


surfaces. 
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Piguce 27 - SEM, TYPLZAL FRACIURE IN SECOND LAY ER ie 
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Piguce 29 - SEM, TYPICAL FRACTURE IN SECOND LAYER, | 208m 
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SES PaeEtEcular Intscest is the“ csontantac, or Srne 
fibers in the layers. All inter-laninarc failur2s occurred 
Game tne epoxy=fiber jolats with n9 discernidle fiber 
breakage. A quantitative analysis, pcesented in Fig. 31, 
indicates that slightly over ten percent of th2 ianer lamina 
epoxy-fiber joints havea broken. These failures are long, 
parallel to the fibers, sontinuous, and appear tS penetrate 
completely. This would indicate anvinabdwiity £5 carry eaay 
load in the transverse jircection while still naintaining 
Semength in the £1lber dicection. This would be sonsistSnt 
With the low transvears2, but very high Wrote lo faticb oliusnew Le 
Samength and failure chacacteristiscs. The transfer of this 
transverse loai to the next lanina would cause large 
meeer=~-laminar shear fo9c2288s, and this could account fer tne 
gross delaminatioos. The observation that all delaminations 
appear to have an ass ciated matrix fracture, or crack, 


further supports this th2ory. 
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2. Microscopic Examiaacion of statics "Yieli" Point 


— =P em «ae == —_ a eee ae oe —_ oP owe = ow 


During static loading of the plates, there was 3 
pronounced "yield" point at approxinately half the ultimate 
road . At this point the panel would experienc2 a decreasad 
load with no increase in jisplacenent, accompaniei by loud 
cracking noises. Figuce2 14 gives plots of thes2 points. A 
plate was loadei just past this "yield" point, and then 
sectioned (5B) as showa in Fig. 13. USiag the erga t 
microscope, Fig. 32 shows an area where all eight lamina are 
delaminated. Ina Fig. 33 thare are just a few ialaminations, 
but both pictures show 2xtensive damage to th2 undersides 


(side opposite th2 load) z2ro degree Layer. 


rigures 34, 35, 113 36 all Show the sana damage to 
the bottom layer along with assorted oth2r delaminations and 
cracks. Figures 34 ani 35 have matrix fractur2s associate2d 
With delaminatiois. Ia the nicrophotographs, the damaged 
outer layer was probably accentuated by the reaction of the 
plastic mounting material with th2 epoxy, as previously 


described. 


Although there w#ere other varied types of damage, 
the major fractire mechanisa was the failure of the zero 
degree lamina on th2 3142 opposite the 1d5ai. Likes 
Gomeecttired “that the “yreld" — point eneountsersi 15 ene 
failure of th2 suter layar in the transverse jlirection and 
the Shifting of this loal t9 tnS innec lananas. Oue to tne 
Secometry of the liy-up ani the Strength charactor stie cs cue 
next two 45-degree laminas could nore efficiently carry the 


load transverse to the zero-degree layer. 


The useful, or sare working, load of th> Sate Would 


be limited by this "yieli" point. 
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Figure 32 - LIGHT MESL9SCOPE, "YESED" BPOGNT EX Ts 
fol PONT Sb yan) OX 
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Figure 33 - LIGHT MNICRISCOPE, "YIELD" POLIT Eee ie een, 
Jo—L LON Soyer 
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Exanination of Initial Dynaaic Damage 


Early in the impact loading history, a "osruise" was 
noted in the imnediat2e indast area. his abnormality could 
mtESt beMdetectsd = using “the Liquids cryse wee fomection 
method. As each test projcressed, the "bruise" could be 
detected by carefully signting across the face of the panel 
where it appeared as a slight bulge. [This danaged area 


never grew larger than th2 immediate area of impact. 


To examine this poenomenon further, a section (1A2) 
was cut from a panel as Shown in Fig. 138. Figure 37 1S in 
an area just outside the 'bruise™ and shows little damages. 
Figure 38 1S in tne edge of the pre-failure danag2 area, ini 
has all plies, with th2 exception of the middle two, 
Srarting "to delaminates. figute 37 Vs Batre ene mp2 ct 


area and delamination is 2xtensive. 


The "bedise”" dslanination and the typical failure 
differ. In the "bruise" area, more lamina ar2 involved 
extensively without a5 nany associated matrix failures in 
adjacent lamina. There ac2 also no Eailures of the inner 


wor lAniNnas»im tne —~ransvVoleSewarescumone 


This damage mechaiism is theorized to pe a form of 
spalling caused by the intensive compressive stress wave 
traveling through the thickness of EhS material on ampact; 


Ref. 10 applies. 
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Pagure 37 = LIGHT NLCRISCOPE, LIGCALIZED DYNAUVIC Data, 
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weer UL CLOSCODIC Exalination Of Baber Panes 


Very few fiber Cailures ware found, and ti2 cause of 


those failures was hard td jieterminso. 


Figure 40 shows a singl2 bcoken fiber in the 45 
degree layer of the Scanning Electron Micross ope sample. 
Figure 41 shows a bundle of broken fibers in th2 zero degrees 


Outer layer of the same specimen. 


The broken fibers in th2se pictures show 45 degree 
fracture surfaces. This would b2 indicative of shear 
E@eeeuce under high stcfaila rates according tomes0ingelom and 
Heiser [Ref. 11). (They sould also just be the 2nds of the 
fibers as manufactured. 


Because 93£ the cacity of beoken fibers ard trove wet 
that fibers are not nocnally continudsus in a Lamina, n9 


Significance was assigned to these fiber fractures. 


It appears the lack of strength of the Laminas in 
the transverse directid1 was the significant factor in 
failure of these slates. Adding 990 degree laminas aignat 
have provided a larj2 tiacraase ina lnpact resistance and 


Static strength. 
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Figure 40° - SEM, A BROKEN @PBER so 270s 
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Seon aonGl REQUIRED FOR FALLURE 


Figure 42 is a plot of the energy required for ultimate 
Pamlure £o0c both dynamic and statis cases. Inly sinply 
Supported plates, using the small penetrator, wec2 plotted. 
A least-squares fit of first and second order wer2 plotted, 
wide Within experimental ascuracy, the enasrgy for failure aa 


both static and dynamic cases was constant. 


Since a constant inpact 2nergy mechanisna was manifest 
and this energy equalled the statics strain energy t9 
failure, the sorrect iynamic nodal seems td be a one 
degree-or-freedom lump2i oarameter spring mass system ([R2f. 
2) oP Tier general,” contLauous, strictucatdynawee wpeoponse 
model was not indicated. However, an extensive series of 
dynamic response studi25 was performed using the dynanic 
response section of SAP [IV [Ref. 5]. Phe parameters of 
impact duration, impact shape, number of natural modes 
Meal zod,;sand Constitutive properties wers all ~ varicde In 
all results th2 2xperimental values showed the actual plate 
Seepeindg ~Significantly stiffer than the "moiel, amie no 
correlation between thea 2xperinental results and the finite 
element dynamic structucal response analysis sould 02 


obtained. 
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BeeeolAkLIC SOLUIIONS 


Wath a constait impact anergy failure, aestatis selutven 
to the plate problen would also provide th2 dynaaic 
Solution . Figuce 8S is a plot of Several “dierorent 


solutions and the average 2xperimental r2sults. 


The exact solution sf Tinoshenk>s and Woinowsky-Kriejer 
Ref. 7] for an orthotropic plate w2area computed aid plotted. 
The design Mmanual's (Ref. 2] exact solution was also plottai 


uSing both sets of daveloo21 constitutive properties. 


The finite element nodel was ausei for the fixed edge 
solution, and the Simply supported solution using the small 
penetrator. It was also used Foc the sinply supported 
solution using th2 larg2 2enetrator. The Jlarge pean2etrator 
was Simulated uSing the 2xperinentally optained 
displacements. A drawinjy was made of the displac2anents with 
the penetrator overlail. Th2 Load waS assuned to be a 
pressure distribution in the sane Shape and proportion to 


the contact area of th2 p2natrator on the displaced plate. 


As can be sean from Fizy. 43, the theoretical results 
agreed well among themsalves. The exact and finite element 
solution for a concentrat2i load wer2 within a percentage 
point; however, the thadretical and 2xperinantal results 


were in significant disagc22ment. 


From the microscopis analysis, tha static theoretical 
analysis, and the dynanisc theoretical analysis it was 
determined that small deflection th2o0ry was unusadle. Thesa 
bates, due tt)  theile snail thickness and relative 


Pectability, were mors appropriately ~Sonstiseede nea 
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deflection theory. 


GChai-Teh Wang Ref. 12] has S@iivyed th]2 Static Srobmoe for 
a large deflection, sinply suppocted, unifornly loaded 
plate. Aithough not directly applicaole t5 the problen 
considered here, his results show th2 san2 general form for 
a uniformly loaded plate as these experimental r2suits for 3 


concentrated load plate. 
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VEE is. CONSEUSTEONS 


Pep ine dynamic and statis failure nodes for simply Supports 


graphite-epoxy plates uni2ac either ceatral concentrated load 


-Or low-energy impact are the same. 


B. Between zero and twenty reet pec second large deflection 


failures of graphite-epoxy, simply supported plat2s occurceid 


at a constant energy. 

Cc. There was a "yield poiat" where danage2 
loading. 

E. There was localized str2ss wave danage 
dynamic contact area. 

F. Linear, small-deflection theory is not 


graphite-epoxy plates loaied to failuce. 
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OCCUELTel In Staci: 


in th2 immediate, 


applicaale sto wenan 








IX. RECOMMENDATIONS 


An exact, larjye deflestion solution should be demeloped 
to analyze thin plate cesdonses to concentrated loads. Once 
this solution is in hand, 1a gooil approxination caa be made 
for low energy, dynamic damage thresholds using energy 
methods. This informatida can then be used for both design 


and operational maintenance. 
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APPENDIX A 
CONSTITUIEVE PRIPERTIES COMPUTER PROGRAM 
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